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Abstract. The region centered in the NGC 6334 1 (N) ra- 
dio continuum source was surveyed in an extension of 6' 
in right ascension and 12' in declination, in the NH3(J,K) 
= (1,1) transition, using the Itapetinga radio telescope. 
The spectra show non-LTE behavior, and gradients of ve- 
locity and line-width were detected along the region. A 
detailed analysis of the spectra showed that the surveyed 
region is composed of at least three overlapped sources 
related to regions that are in different stages of star for- 
mation: NGC 6334 1, associated with an already known 
molecular bipolar outflow, NGC 6334 1(N)w, the brightest 
ammonia source, coincidental with the continuum source 
NGC 6334 1(N), and NGC 6334 1(N)e, weaker, more ex- 
tended and probably less evolved than the others. The 
physical parameters of the last two sources were calcu- 
lated in non-LTE conditions, assuming that their spectra 
are the superposition of the narrow line spectra produced 
by small dense clumps. The H2 density, NH3 column den- 
sity, kinetic temperature, diameter and mass of the clumps 
were found to be very similar in the two regions, but the 
density of clumps is lower in the probably less evolved 
source NGC 6334 1(N)e. Differences between the physical 
parameters derived assuming LTE and non-LTE condi- 
tions are also discussed in this work. 

Key words: stars; formation - ISM: individual (NGC 
6334) - ISM: abundances, clouds, H 11 regions, molecules 
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1. Introduction 

The NGC 6334 region, located at a distance of 1.74 
kpc from the Sun (Neckel |l978|) , is a c omple x of ra- 
dio an d inf rared sources (Sh aver fc Goss 1970 ; Haynes 
et al. 



1986) but it is below the detection limit at wavelengths 



1975; McBreen et al. 1979). One of its compo- 



nents, NGC 6334 1(N), is also the strongest source in th e 
NH3(J,K) = (1,1) transition (e.g., Schwartz et al. |1978D . 
It is also a strong co ntinuu m sourc e at 4 00, 1000 and 1300 
fim (Cheung et al. 1978 Gezari 1982| ; Loughran et al. 



smaller than 137 iim. It has an II2O maser source, CHaOH 



masers (Mo ran fc Rodriguez 1980 Kogan fc Slysh 1998; 
Walsh et al. 199S), six reddened sources at the JHK bands 



(Tapia et al. 1996), a bipolar molecular outflow in Si O an d 
eight knots of H2 emission (Megeath & Tieftrunk 1999). 
However, no H 11 region associated to this source was de- 
tected. 

The source NGC6334 I, located approximately 2' south 
of NGC6334I(N), h as H 2O, OH and CH^H masers 



(Moran fc Rodr iguez |1980[ Gaumefc Mutel |1987| ; Forster 
fc Caswell I989[ Menten fc Batrla 198£ ) , an ultra-compact 



Hii region (Rodriguez et al. 1982 ; DePree et al. 1995 ) and 
a c luster of stars detected at the JHK bands (Tapia et 
al. 1996). It was studied with high angular resolution in 
the NH3(J,K) = (1,1) transition by Jackson et al. (|1988|) 



(hereafter, JHH88), who interpreted the observed struc- 
ture as a circumstellar molecular disk rotating around a 30 
Mq O star. On the other hand, CO and H2 observations 
suggest that the structure arises from a molecular bipola r 
outflow (Bachiller & Cernicharo 1990; Persi et al. 1996). 



Due to these characteristics, it is probable that this re- 
gion is in a more advanced stage of stellar formation than 
NGC 6334 1(N). 

Both regions have been studied in several transitions 
of i^CO, "CO and CS (Kraemer & Jackson |1999D , the 
meta-stable transitions of the ammonia molecule, includ- 
ing NH3(J,K) = (1,1), ( 2,2), ( 3,3) and (6,6) [e. g., Sch wartz 
et al. |1 9 78| ; Fors ter et al. [1987| Vilas-B oas et al. |1988t Jack- 
son et al. 1988| ; Kuiper et al. 1995; Kraemer fc Jackson 
|1995 , 199£ I, and the non-meta-stable transition NH3(J,K) 
— (2,1) (Kuiper et al. 1995). The physical parameters of 
these regions were always derived from the NH3(J,K) = 
(1,1) and (2,2) observations, assuming LTE conditions. In 
fact, the intensity anomalies in the nuclear quadrupole hy- 
perfine structure of the NH3(J,K) ~ (1,1) spectrum are 
usually neglected, although they are observed in many 
warm molecular clouds (Stutzki et al. 1982; Batrla et al. 



1983t Stutzki et al. |l984D 
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Matsakis et al. ( 1977 ) interpreted the anomalous spec- 
trum as the consequence of non-thermal population in 
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the hyperfine states, induced by selective trapping in the 
hyperfine transitions of NH3(J,K) = (2,1) (1,1). This 
effect is relevant only when the width of the hyperfine 
lines lie between 0.3 and 0.6 km s~^. Due to this restric- 
tion, they assumed that the molecular cloud is formed by 
clumps, each one producing a narrow line spectrum, so 
that the observed spectrum would be the superposition 
of individual clump spectra. Based on this model, Stutzki 
& Winnewisser ( 1985D (hereafter, SW85) elaborated a nu- 
merical algorithm to calculate the brightness temperature 
ratios between the hyperfine satellites and the main line 
in the NH3(J,K) = (1,1) spectrum, and also the bright- 
ness temperature ratio between the (2,2) to the (1,1) main 
lines, as a function of the NH3 column density, H2 density 
and kinetic temperature of each clump. 

Gaume et al. ( 1996| ) questioned this model, based on 
observations of the NH3(J,K) = (1,1) absorption spectra 
towards the continuum source DR21. Their argument was 
based on an anti-correlation in the degree of LTE depar- 
ture of the inner and outer hyperfine components in differ- 
ent regions of the cloud. They proposed instead, inflows 
and outflows of matter to explain the observations, but 
did not present any quantitative calculation to support 
their claim. 

In this work, we used our high quality observations 
of the NGC 6334 region in the NH3(J,K) = (1,1) transi- 
tion, which allowed the precise determination of the inten- 
sity anomalies in the hyperfine structure, to obtain with 
good reliability its physical parameters assuming non-LTE 
conditions. The data allowed us also to separate in the 
spectra the contribution of at least three different sources: 
NGC 6334 1, NGC 6334 1(N)w and NGC 6334 1(N)e. 

2. Observations 

The observations were made during July 1996 using the 
13.7 m radome enclosed radio telescope of the Itapetinga 
Radio Observatory]]. The receiver front-end consisted of a 
circularly polarized corrugated horn connected to a cooled 
HEMPT. An acusto-optical spectrometer was used in the 
back-end, with spectral resolution and total bandwidth of 
70 kHz and 41 MHz respectively. The total system temper- 
ature was about 200 K. The angular resolution of the radio 
telescope at the frequency of NH3(J,K) — (1,1) transition 
was 4.2'. The observations were made using the ON-OFF 
total-power technique, switching between positions every 
20 s with amplitude of 20' in azimuth, enough to guar- 
antee that there was no source in the OFF position. The 
effective integration time of each observation was 210 s. 
A 15 K noise source and a room temperature load were 
used in the calibration to obtain the gain and to correc t 
for atmospheric attenuation (Abraham & Kokubun 1992|) . 
The continuum point source Virgo A was used as a pri- 
mary calibrator. The NGC 6334 region was mapped in 

^ Operated by CRAAE, Centro de Radio Astronomia e 
Aplicagoes Espaciais, Sao Paulo, Brazil 



the NH3(J,K) = (1,1) transition, covering 6' in right as- 
cension and 12' in declination, with a spacing of 2'. The 
central position of the map, with equatorial coordinates 
a(1950) = 17'^17'"32" and 5(1950) = -35°42', was ob- 
served several times during the mapping period and was 
used as secondary calibrator. This position was also ob- 
served in the NH3(J,K) = (2,2) transition. 

A polynomial baseline was subtracted from the obser- 
vations and the five Gaussian functions were fitted to the 
hyperfine NH3(J,K) — (1,1) spectra. The velocity separa- 
tion between the Gaussians was fixed and the line-width 
of the hyperfine satellites constrained to be equal to the 
main line. This last assumption is well justified, since the 
optical depth of the main line is small (t ~ 1.4). 

The results from the observations are presented in Ta- 
ble 1. The first column shows the position in the map 
relative to the central coordinate, columns 2-6 exhibit the 
antenna temperature of each transition corrected for at- 
mospheric attenuation, columns 7 and 8 the velocity and 
line- width of the main line and the last column the num- 
ber of observations made in each position. The NH3(J,K) 
= (1,1) spectra of some points of special interest, as well 
as the NH3(J,K) = (2,2) spectrum of the central position 
are shown in Figure 1. 

The NH3(J,K) = (1,1) spectra in all the observed po- 
sitions exhibit systematic differences between the line in- 
tensities in each pair of hyperfine transitions (F=l— ^O, 
F=0^1) and (F=1^2, F=2^1), which are, at least in 
the strongest points, much larger than the rms measured 
at the baseline. These intensity differences can be inter- 
preted as the signature of non-LTE conditions. However, 
we should point out that this conclusion cannot be ob- 
tained by looking only at the residuals of the line fittings, 
which present larger rms than the baseline, both under 
LTE and non-LTE assumption. This is expected an ex- 
pected result, since the superposition of sources and the 
existence of velocity gradients can make the line shapes 
differ from Gaussian functions. 

The antenna temperature distribution and gradients 
of velocity and line- width are presented in Figure 2. They 
agree with t he results of Forster et al. (1987) and Kuiper 
et al. (1995), who mapped with high spatial resolution a 
small region (about 100" in diameter) around the central 
position in our map in the NH3(J,K) = (3,3) transition. 
However, the distribution of velocities differs from that 
found by Schwartz et al. (197S), since their map do not 



show an increase in the velocity's modulus towards neg- 
ative declinations, possibly due to the low signal-to-noise 
ratio in their spectra. 

3. Results 

3.1. Identification of the sources NGC 6334 I, 
NGC 6334 HN)w and NGC 6334 I(N)e 

Although the map in Figure 2 seems to be consistent with 
the presence of a single source centered in NGC 6334 1(N), 
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Table 1. Parameters obtained from the Gaussian fitting for each position in the map. 



(Aa, Aa)" 


Ta(F=1^0) 


Ta(F=1^2) 


Ta(AF=0) 


Ta(F=2^1) 


Ta(F=0^1) 


vlsr ^ 


AvpwHM 


#° 


(arcmin) 


(K) 


(K) 


(K) 


(K) 


(K) 


(km s"^) 


(km s"-') 




(-2,-2) 


0.17 ± 0.04 


0.27 ± 0.04 


0.67 ± 0.04 


0.24 ± 0.04 


0.26 ± 0.04 


-5.09 ± 0.12 


6.00 ± 0.13 


4 


(-2,0) 


0.15 ± 0.03 


0.24 ± 0.03 


0.54 ± 0.03 


0.21 ± 0.03 


0.20 ± 0.03 


-4.40 ± 0.12 


4.28 ± 0.12 


4 


(-2,2) 


0.09 ± 0.03 


0.13 ± 0.03 


0.43 ± 0.03 


0.15 ± 0.03 


0.14 ± 0.03 


-3.82 ± 0.12 


3.06 ± 0.13 


3 


(0,-4) 


0.15 ± 0.04 


0.23 ± 0.04 


0.39 ± 0.04 


0.12 ± 0.04 


0.11 ± 0.04 


-4.66 ± 0.12 


6.29 ± 0.12 


4 


(0,-2) 


0.30 ± 0.03 


0.43 ± 0.03 


0.86 ± 0.03 


0.33 ± 0.03 


0.32 ± 0.03 


-4.98 ± 0.12 


5.57 ± 0.12 


2 


(0,0) 


0.41 ± 0.01 


0.62 ± 0.01 


1.41 ± 0.01 


0.53 ± 0.01 


0.52 ± 0.01 


-4.66 ± 0.12 


4.26 ± 0.12 


19 


(0,2) 


0.30 ± 0.04 


0.51 ± 0.04 


1.05 ± 0.04 


0.42 ± 0.04 


0.43 ± 0.04 


-4.21 ± 0.12 


4.06 ± 0.12 


2 


(0,4) 


0.15 ± 0.03 


0.24 ± 0.03 


0.62 ± 0.03 


0.17 ± 0.03 


0.29 ± 0.03 


-3.36 ± 0.12 


3.11 ± 0.12 


4 


(2,-4) 


0.10 ± 0.04 


0.22 ± 0.04 


0.44 ± 0.04 


0.21 ± 0.04 


0.21 ± 0.04 


-5.98 ± 0.12 


5.37 ± 0.12 


2 


(2,-2) 


0.27 ± 0.04 


0.39 ± 0.04 


0.86 ± 0.04 


0.34 ± 0.04 


0.34 ± 0.04 


-5.29 ± 0.12 


4.57 ± 0.12 


2 


(2,0) 


0.38 ± 0.02 


0.54 ± 0.02 


1.28 ± 0.02 


0.47 ± 0.02 


0.48 ± 0.02 


-4.74 ± 0.12 


4.09 ± 0.12 


6 


(2,2) 


0.34 ± 0.06 


0.32 ± 0.06 


1.06 ± 0.06 


0.47 ± 0.06 


0.43 ± 0.06 


-4.43 ± 0.12 


4.21 ± 0.13 


1 


(2,4) 


0.16 ± 0.03 


0.28 ± 0.03 


0.59 ± 0.03 


0.25 ± 0.03 


0.23 ± 0.03 


-3.71 ± 0.12 


3.53 ± 0.13 


3 


(2,6) 


0.13 ± 0.03 


0.18 ± 0.03 


0.38 ± 0.03 


0.19 ± 0.03 


0.16 ± 0.03 


-3.07 ± 0.12 


2.24 ± 0.13 


4 


(2,8) 


0.10 ± 0.04 


0.10 ± 0.04 


0.25 ± 0.04 


0.07 ± 0.04 


0.10 ± 0.04 


-3.42 ± 0.12 


4.10 ± 0.12 


4 


(4,-2) 


0.04 ± 0.03 


0.07 ± 0.03 


0.19 ± 0.03 


0.05 ± 0.03 


0.08 ± 0.03 


-4.85 ± 0.13 


3.95 ± 0.14 


4 


(4,0) 


0.15 ± 0.03 


0.18 ± 0.03 


0.47 ± 0.03 


0.14 ± 0.03 


0.12 ± 0.03 


-4.75 ± 0.12 


4.09 ± 0.13 


8 


(4,2) 


0.10 ± 0.03 


0.14 ± 0.03 


0.35 ± 0.03 


0.12 ± 0.03 


0.13 ± 0.03 


-4.26 ± 0.12 


4.23 ± 0.12 
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Offset in 


relation to the 


equatorial coordinates a(1950) — 


17''17"32° and 5(1950) = -35° 


42'. 







^ Velocity of the main line. 
° Number of observations made in each position. 



high resolution observations (JHH88; Kraemer & Jackson 



199E) showed the existence of at least two ammonia fea- 



tures, with angular sizes of the order and smaller than 
1', centered at the positions of the continuum sources 
NGC6334I(N) and NGC6334I respectively. Taking into 
account the half-power beam width of the Itapetinga radio 
telescope, we find that the expected contribution of these 
two sources at the position (Aa, A6) = (2', 6') [hereafter, 
just (2', 6')] should fall below the detection Hmit given by 
the rms of the observations. Since this is not the case, 
as can be seen in Figure 1(d), we conclude that another 
source is present in the region. We will describe, in the 
following paragraphs, the procedure used to separate the 
contribution of each component. 



NGC 63341, when observed with high resolution 
(JHH88), exihibts two main structures, with velocities 
centered in -6.6 and -9.4 km s~^ and widths of 2.0 and 3.5 
km s~^ respectively. The maximum flux densities of these 
structures are 0.42 and 0.48 Jy beam~^ respectively which, 
once diluted into the 4.2' Itapetinga radio telescope beam, 
correspond to antenna temperatures of 0.20 and 0.23 K re- 
spectively. We calculated the contribution of this source 
to the NH3(J,K) = (1,1) spectra in our map as the con- 
volution of the antenna beam and two point sources, with 
the characteristics mentioned above. The intensity of the 
hyperfine satellite lines was calculated assuming LTE in 
an optically thin source. Although the antenna tempera- 
ture of NGC 6334 1 is small compared to the intensity of 
NGC6334I(N), its effect on the line-width is significant. 



and its subtraction decreased the observed gradient along 
the map. 

Since the high resolution observations of Kraemer & 
Jackson ( |1999D have revealed NGC 6334 1(N) as an ex- 
tended source, although smaller than the Itapetinga radio 
telescope beam, it was difficult to subtract its contribu- 
tion to expose the third and much weaker source. For this 
reason we used the spectrum at position (2', 6'), where 
we did not expect any contribution from the strongest 
source, as representative of the spectrum of the weak- 
est. We labeled these sources as NGC 6334 1(N)w and 
NGC6334I(N)e [hereafter, I(N)w and I(N)e] respectively. 
We assumed no velocity gradients and subtracted a frac- 
tion of the (2', 6') spectrum from all the others. This frac- 
tion was chosen in such a way that the resultanting spec- 
tra had the same departure from LTE as found at position 
(0',0'). Moreover, since there is a velocity gradient along 
the region, the subtracted fraction was also limited in or- 
der to avoid fake absorptions in the resulting spectra. 

In Figure 3, we show the superimposed antenna tem- 
perature maps of the three regions, NGC 6334 1, 1(N)w and 
I(N)e. We can see that I(N)w is aproximately 2.5 and 6.5 
times brighter than I(N)e and NGC 6334 1, respectively. 
The source I(N)w turned out to be smaller than the beam 
size, as expected from the observations of Kraemer & Jack- 
son (1999). On the other hand, the I(N)e region looks more 
elongated in declination. 

The subtraction procedures affected the magnitude of 
the gradients of velocity and line-width, but not their di- 
rections. In fact, the gradients became softer than before. 
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Fig. 1. Spectra from the central (Aa, AJ) = (0',0') position for the transitions (a) NH3(J,K) = (1,1) and (b) (2,2). 
The NH3(J,K) = (1,1) spectrum from the (c) (0',-2') and (d) (2', 6') positions are also shown. The thin lines represent 
the Gaussian fittings, while the lines, below the observed spectra, are the residuals. 



showing that the other sources were responsible in some 
amount for the variations of these properties along the 
region I(N)w. As it can be seen in Figure 4, the veloc- 
ity gradient is aproximately perpendicular to the Galactic 
plane. According to the model adopted here, the gradient 
in line-width is related to gradient in turbulence, which 
turned out to be larger in the direction of the more evolved 
sources. 

3.2. The physical parameters of I(N)w and I(N)e regions 
in non-LTE conditions 

We used the antenna temperature of the hyperfine satel- 
lites relative to the main line in the NH3(J,K) = (1,1) 
spectrum and the ratio between the antenna temperatures 
of the main lines of NH3(J,K) = (2,2) and (1,1), shown 
in Table 2, to calculate the relevant physical parameters 
for the I(N)w and I(N)c regions. We applied the model 
of SW85 to obtain the NH3 column density {N^Hs), the 
kinetic temperature (Tk) and the H2 density (riHs) of the 
clumps in the two regions. In region I(N)w we used the 
spectrum from the (0',0') position, which corresponds to 



the highest intensity, and in I(N)e we used the (2', 6') posi- 
tion, where we believe there is not important contribution 

from I(N)w. 

To determine the physical parameters of the clumps 
we looked into the (TVnHs, J^Hai ^k) space for the point 
at which, according to the model of SW85, the curves 
representing the observed line intensity ratios intercept 
themselves. We had to extrapolate these curves up to H2 
densities of 10^ cm~^ to obtain a good crossing point, 
which occured at ~ 10^ ''' cm~^. 

In Figure 5, we show the curves of constant ratio be- 
tween the hyperfine satellites and the main line in the 
NH3(J,K) = (1,1) spectrum, as well as the ratio be- 
tween the main line intensities of the NH3(J,K) = (2,2) 
and (1,1) transitions for a kinetic temperature of 24 K, 
which corresponds to the best fit in both regions. In Ta- 
ble 3, we present the values of the physical parameters 
extracted from this figure: the brightness temperature for 
the NH3(J,K) = (1,1) and (2,2) transitions [Tb(1,1) and 
Tb(2,2), respectively], the kinetic temperature, the NH3 
column density and the H2 density. We can see that the 
physical conditions of the clumps in both regions are prat- 
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Fig. 2. Maps of (a) antenna temperature (in units of K), (b) velocity (in units of km s ^) and (c) line-width (in units 
of km s-i) for the NGC 6334 region. The (0',0') position has coordinates a(1950) = 17''17™32'* and 5(1950) = -35°42'. 
The black dots mark the observed positions. 



Table 2. Ratio between the observed antenna temper- 
atures of each hyperfine satellite in relation to the (1,1) 
main line [T™] , besides the ratio between the antenna tem- 
perature of the (2,2) and (1,1) main lines. 

NGC 6334I(N)w NGC 6334 1(N)e 



Ta(f=i^o) / rr 


0.295 ± 0.008 


0.34 ± 0.09 


Ta(f=i^2) / rr 


0.409 ± 0.008 


0.47 ± 0.09 


Ta(F=2^1) / Tr 


0.371 ± 0.008 


0.50 ± 0.09 


Ta(F=0^1) / 


0.364 ± 0.008 


0.42 ± 0.09 


Ta(2,2) / rr 


0.67 ± 0.04 





tra, their dispersion velocities must be small, implying 
that they should be gravitationally stable. Therefore, their 
individual masses (Md) cannot be larger than the Jeans 
mass (Mj) and their diameters (Ld) must be smaller or 
equal than the Jeans length {Lj). Arguments in favor of 
t hese hypotheses have been presented by Curry & McKee 
( ^OOOl) based on numerical simulations. Assuming spheri- 
cal clumps with a ratio of helium to molecular hydrogen 
of 0.25, we calculated their diameter and mass, using the 
following expressions: 



Lr] — i.T — 



TT k Tk 



(1) 



ically the same. It is important to stress that the val- 
ues listed in Table 3 would not change substantially if 
we use the original spectra, including the contribution of 
the three sources, since I(N)w is much brighter than the 
others. On the other hand, the dispersion of the contours 
around the intersection point in Figure 5(a) decreased ap- 
preciably after the separation of the three components, 
resulting in a better accuracy in the determination of the 
physical parameters of the region. 

3. 3. Physical parameters of the individual clumps 

As it was previously mentioned, to explain the hyper- 
fine structure intensity anomalies in the NH3(J,K) = (1,1) 
transition by the superposition of individual clump spec- 



Mcl = A/, J = -mn riHa -^ci 



(2) 



where k is the Boltzmann's constant, mn the mass of hy- 
drogen atom and G the gravitational constant. 

The beam filling factor (/jk) can be calculated from 
the brightness temperature, Tb{J,K), derived from the 
model and T^^{J,K)^ obtained from the observations, 
which is proportional to the antenna temperature: 



/jK — 



T^'%J,K) ^ 2 c' A Ta{J,K) 
Tb(J,X) TTku'elTB{J,K) 



(3) 



where 9^ is the half-power beam width of the telescope 
and A is the ratio between flux density and antenna tem- 
perature of a continuum point source. In this work, we 
used the radio galaxy Virgo A as the calibrating source. 
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Fig. 3. The superposition of the contribution of the three 
sources: NGC 63341 (dotted hnes), I(N)w (sohd hnes) and 
I(N)e (dashed hnes). The contours are given in antenna 
temperature, starting at 0.1 and ending at 1.3 K (with 
spacing of 0.1 K). The gray points mark the observed po- 
sitions. The rectangle shows the area that was observed 



by Tapia et al. 199(: in the JHK bands. The black stars 



are the sources which are detected by these authors. 

Using the definition of beam filling factor, we can de- 
termine the number Kd of clumps inside the beam's solid 
angle: 

Kd=hK(^y (4) 



where D is the distance from the cloud. 

To obtain the density of clumps (Cci), we can consider 
two different situations. The first one corresponds to the 
clumps spread over a solid angle larger than the beam. In 
this case, the density of clumps is given by: 

4 K,i 

Ccl = - 




♦4.0 43.0 *2.0 +1.0 0.0 -1.0 -2.0 
Right Ascension Offset (arcmin) 




+4.0 +3.0 +2.0 +1.0 0.0 -1.0 -2.0 
Right Ascension Offset (arcmin) 

Fig. 4. The gradients of (a) velocity and (b) line-width 
for region I(N)w. The black points mark the observed po- 
sitions. The countour values are given in units of km s^^. 



where 9^ and Op are, respectively, the extension of the re- 
gion in right ascension and declination. In the second one, 
if the source is smaller than the antenna beam, we can 
only estimate either an upper limit, considering that the 
clumps are grouped side by side and without superposi- 
tion, or a lower limit, supposing they are spread over the 
whole beam: 



6 K, 



cl 



(5) ^{DOi 



< Ccl < -- 



(6) 
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NGC6334 l(N)w - T,^ =24 K 




5.0 6.0 7.0 

log[n(H2) (cm-3)] 

NGC6334 l(N)e - 1^=24 K 



— 15.8 
E 

CO 

CM 

E 15.4 

X 

S 

1^15.0 

n: 



(b) 



14.6 



a) 
o 



14.2 



4.0 



5.0 6.0 

log[n(H2) (cm-3)] 



7.0 



Fig. 5. Ratio between the intensity of the sateUites and 
the main hne of the NH3(J,K) = (1,1) spectrum, rep- 
resented by: dashed, dashed dotted, dashed dotted dot- 
ted and dotted hnes (corresponding to F=l— >0, F=l— >2, 
F=2^1 and F=0^1 transitions, respectively) for a ki- 
netic temperature of 24 K. The upper and lower plots 
[(a) and (b) respectively] are associated to the I(N)w and 
I(N)e in this order. Ratio between the main lines intensity 
of the NH3(J,K) — (2,2) and (1,1) transitions is shown 
by the solid lines. The weaker lines are the uncertainties 
associated to the fitting. The locus where these curves in- 
tersect each other represents the physical parameters of 
the region. 

Once the density of clumps and their diameters are 
known, we can calculate their collision time (tcoii): 



coll 



CciAi;j^^ 



(7) 



FWHM 



We can also calculate the mean H2 density [fiYi^, 
which is the density of an homogeneous cloud with mass 
equal to the total mass of the clumps filling the telescope 



Table 3. The physical parameters determined from the 
Figure 5 for I(N)w and I(N)e regions. 





NGC6334I(N)w 


NGC6334I(N)e 


rB(i,i) (K) 


12.6 ± 0.4 


15.5 ± 1.0 


Tb(2,2) (K) 


8.4 ± 0.6 


11.6 ± 1.1 




24.0 ± 0.8 


24.0 ± 1.0 


AfNH3 (10'"^ cm-2) 


2.3 ± 0.8 


2.4 ± 0.8 


nn^ (10^ cm-3)->b 


1.8 ± 0.4 


1.4 ± 0.7 



For each clump. 

According to SW85, the H2 density derived from their 
model must be divided by 1.75 in order to obtain the true 
value for this parameter, shown in this table. The reason 
for this procedure is that only NH3 - He collisions were 
considered in the numerical model of SW85 as excitational 
mechanism of the ammonia molecule. 



beam: 



cl 



(8) 



3 mn 

The ammonia abundance is calculated from the follow- 
ing expression: 



X 



NH3 



-^cl»^H2 



(9) 



where it was assumed that there is not superposition be- 
tween the clumps, so that the line of sight extension is 
equal to the diameter of each clump (SW85). 

The physical parameters of the clumps in regions 
I(N)w and I(N)e and their respective uncertainties are 
listed in Table 4. The quoted errors for the mass and di- 
ameter of the clumps represent the uncertainties in the 
kinetic temperature and H2 density. 

The clumps in both regions have similar diameters 
(0.007 pc), masses (0.2 Mq) and ammonia abundances 
(7xl0~^). They also seem to be gravitationally bounded, 
since the total virial mass is about a factor of ten lower 
than the total mass of the clumps which are contained 
in the beam's sohd angle [7600 and 2100 for I(N)w 
and I(N)e respectively]. The 36000 clumps of I(N)w are 
spread in such a way that their density is between 6000 
and 33000 clumps per pc~^, much larger than the 720 
clumps per pc~'^ in the I(N)e region. 

The beam filling factor for the I(N)w region is larger 
than expected from the NH3(J,K) — (3,3) data presented 
in Kraemer & Jackson ( 1999D (around 0.10 for a beam size 
of 4'). The reason is that we had assumed no superposition 
between the clumps along the line of sight. The physical 
superposition required by the smaller filling factor is small 
enough to guarantee very little velocity superposition, in 
which case the spectra of the superposed clumps can be 
simply added. 
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Table 4. The physical parameters for I(N)w and I(N)e re- 
gions calculated in non-LTE conditions. The uncertainties 
are given in parenthesis. 





NGC6334I(N)w 


NGC 6334I(N)c 


id (10-^ pc) 


6.7(0.8) 


7.5(1.7) 


Moi (Mq) 


0.21(0.03) 


0.24(0.06) 


/jK 


0.33(0.03) 


0.10(0.01) 


Kci (10^) 


36(10) 


9(4) 


Cci (10=* pc-^) 


6(2) < Col < 33(11) 


0.72(0.36) 


tcoll (10° years) 


4(2) < tooii < 22(9) 


230(155) 


(10* cm"^) 


1.7(0.6) < fiHa < 9.5(3.4) 


0.23(0.13) 


XnH3 (10-^) 


6.2(2.3) 


7.2(3.0) 


A/totClO^MQ)" 


7.6(2.3) 


2.1(1.2) 


A/vir(lO''AfQ)'' 


0.67(0.06) 


0.15(0.02) 



Total mass of whole clumps which are contained in the telescope 
beam. 

Mass obtained from the Virial Theorem. 



4. Discussion 

4-.1. Ammonia abundance for the regions I(N)w and 
I(N)e 

The study of ammonia abundance is important because it 
provides observational constrains for the chemical models 
that try to explain how this molecule is formed in the in- 
terstellar medium, and it may be used as an age indicator 
of the molecular cloud (e.g., Suzuki et al. 1992 ). Although 
the chemical models calculate the evolution of the ammo- 
nia abundance as a function of time with a limited confi- 
dence, due mainly to the uncertainties concerning the for- 
mation of this molecule in the interstellar medium, they 
agree that the time necessary to reach the equilibrium 
abundance is approximately 10^ years. If depletion effects 
are not consider in the model, the equilibrium abundance 
for the ammonia molecule is 10~^. If there is depletion, 
this value might be 10^^ (Bergin et al. 1995| ) or 10~^ 



(Hasegawa & Herbst 1993 ), even though in this last case 
the abundance reaches equilibrium after 10* years. 

Comparing the ammonia abundances listed in Table 4 
with the equilibrium abundance, we conclude that either 
the regions I(N)w and I(N)e have an age lower than 10^ 
years, or that ammonia depletion is important in both 
sources. This result was also obtained by Kuiper et al. 
(|l99|), assuming LTE conditions. Another relevant aspect 
that must be taken into account is the dependence of the 
equilibrium abundance on the nitrogen abundance in the 
interstellar medium, and therefore the gradient of nitro- 
gen with galactocentric distance must be also considered 
(Vilas-Boas & Abraham pOOO]). 



4-. 2. Comparison with LTE models 

In this sub-section, we will derive, under assumptions of 
LTE, the physical parameters of the source I(N)w and 



compare them with the parameters derived in section 3 
under non-LTE conditions. We calculate these parameters 
at the (0',0') position because it is the only point in which 
the NIl3(J,K) — (2,2) transition was observed. In the LTE 
calculation we considered a three-level system, in order to 
obtain the rotational and kinetic temperatures (Tr and Tk 



respectively) through the equations (Ho & Townes 1983 
Walmsley & Ungerechts 119831): 



-41.5 X 



In 



9 T22 Aw22 



FWHM 



20 Tn 



A?; 11 



l + (TK/41.7ii:)ln 



1 



C(22^21) 



C(22^11) 



(10) 



(11) 



where rn and T22 are respectively the optical depths in 
the transitions NH3(J,K) = (1,1) and (2,2), C(22 21) 
and C(22 1 1) ar e the collision rate coefficients given by 
Danby et al. ( |l988D . 

To determine the optical depth, we assumed that the 
cloud has a spherical shape, that the excitation tempera- 
ture (Tex) is the same in all hyperfi ne lin es of the NH3(J,K) 
= (1,1) spectrum (Ho & Townes 1983 ). Thus, using the 
ratio between the (1,1) main line and its inner hyperfine 
satellites, rn can be calculated from: 



Tl-(1,1) l-e(Tfi) 



?a(1,1) 



1 



(12) 



where r™ and t^-^ are, respectively, the optical depth for 
the main line and the inner satellites, with rfj^ = 0.28 Tf{ 
(Ho & Townes |l98§) and e(r) = 2t-^[1 - e-^(l -I- r)] 
(SW85). To obtain the optical depth for the NH3(J,K) 
— (2,2) transition, it is necessary to substitute in (12) 
T^(l, 1) by Ta(2, 2), the antenna temperature of the (2,2) 
main line. 

As we have less equations than unknown parameters, 
it is necessary to make some assumptions to solve the sys- 
tem. There are at least two clear alternatives: we may 
assume that either the excitation and rotational temper- 
atures have the same value, or that the source fills com- 
pletely the antenna beam. Here, we adopted the first al- 
ternative since high resolution observations indicate fill- 
ing factors smaller than one for these sources (Kuiper et 



al. 


1995 


1999 


)• 



We present in Table 5 the physical parameters ob- 
tained under LTE and non-LTE conditions. The NH3 col- 
umn density and the ammonia abundance are not sen- 
sitive to LTE departures. The largest discrepancies are 
found in the kinetic and rotational temperatures. In fact, 
the kinetic temperature estimated under LTE conditions 
is about twice the value obtained in non-LTE. The opti- 
cal depth of I(N)w, in both transitions, is slightly smaller 
than what it was found in previous studies (e.g., Forste r 
et al. [I987t Vilas-Boas et al. |l988| ; Kuiper et al. [19951 ). 
We believe that the differences can be attributed to the 
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Table 5. Comparison between the physical parameters 5. Conclusions 



for I(N)w 


calculated in LTE and 


non-LTE conditions. 




LTE 


non-LTE 


(K) 


44 ± 8 


24.0 ± 0.8 


Tr (K) 


28 ± 3 


18.7 ± 0.3 


Til 


1.42 ± 0.01 




r22 


0.71 ± 0.08 




fll 


0.88 ± 0.10 


0.33 ± 0.03 


f22 


0.88 ± 0.10 


0.33 ± 0.04 


fiH2 (10* cm-^) 12 ± 5 1.7 


± 0.6 < < 9.5 ± 3.4 


A'NHg (10^^ cm^^) 2.5 ± 0.1 


2.3 ± 0.8 


XnH3 (10 


3.3 ± 1.7 


6.2 ± 2.3 



lower signal-to-noise ratio in previous observations. The 
beam filling factor is larger under LTE conditions. This 
was already expected because under non-LTE conditions 
the gas is confined in clumps and not spread in a homo- 
geneous cloud. 

4-.3. Evolutionary state of the regions 

We will try to determine now the evolutionary state of 
the three regions identified in our ammonia observations. 
A possible way to explore this question is comparing the 
distribution of indicators of early stages of star forma- 
tion. In Figure 3, we show the distribution of infrared 



sources, detected in the JHK bands by Tapia et al. ( |1996D 
in a small region of the mapped area. We can see a de- 
crease in the number of infrared sources towards the re- 
gion I(N)w and also the existence of a cluster of sources 
around NGC 6334 1, indicating that it is more evolved 
than I(N)w. Other observational results favorable to this 
scenario are the detection of an UCHii region and the 
identification of a molecular bipolar outflow associated to 



this source (Rodriguez et al. 1982; Bachiller & Cernicharo 
199C). There are no infrared observations towards I(N)e, 
to indicate its evolutionary stage. 

Our observations, together with the non-LTE model 
(SW85), show that the physical properties of the clumps 
in I(N)w and I(N)e are similar. However, the volume den- 
sity of clumps in I(N)w (> 6000 pc~'^) is at least a hun- 
dred times higher than in I(N)e, resulting in a collision 
rate which is hundred times higher in I(N)w than I(N)e. 
Thus, assuming that collisions result in coagulation of the 
fragments, which finally collapse to form stars after sev- 



eral collision times (SW85; Scalo & Pumphrey 1982), we 
expect that I(N)e is in a earlier star formation stage than 
I(N)w. 

The evolution scenario proposed by SW85 and the 
results presented in this work indicate that the regions 
I(N)w and I(N)e could form in the future star clusters, 
similar to those detected in NGC 6334 1. 



The NH3(J,K) = (1,1) spectra from the NGC 6334 region 
exhibit hyperfine structure intensity anomalies, suggest- 
ing departures from LTE conditions. We have used the 
numerical model o f SW 85, based on the qualitative model 
of Matsakis et al. (1977) to calculate the physical parame- 
ters of the ammonia clouds. This model takes into account 
the non-LTE effects through selective trapping in the far 
IR transition NH3(J,K) = (2,1) — > (1,1) and hyperfine se- 
lective collisional excitation. Its main limitation is that it 
does not include the effects of an IR continuum at the 
NH3(J,K) = (2,1) (1,1) transition frequency. As it was 
stressed by SW85, its influence on the anomalies is diffi- 
cult to predict, but they argued that if the IR radiation 
field had significant influence, variations in their strength 
with the distance from the exciting source would be seen, 
which is not supported by observations, at least in the 
regions S106 and W48. 

The existence of clumps with very small dispersion ve- 
locity was questioned by Gaume et al. ( 1996 ), based on the 
NII3 absorption spectra in the direction of the continuum 
source DR21. They proposed, instead, inflows and out- 
flows of matter in different regions of the cloud. However, 
the mechanism producing the non-LTE level population 
is the same as in Matsakis et al. ( 1977 ), that is, the selec- 
tive trapping of photons in the NH3(J,K) = (2,1) (1,1) 
transition. Since this model was not developed quantita- 
tively, it is not clear if it will work to explain the observed 
anomalies. 

Although our observations do not have high angular 
resolution, our high quality spectra allow the identifica- 
tion of three distinct sources. The first one is NGC 6334 1, 
a compact source mapped in NH3(J,K) = (1,1) transi- 
tion by JHH88, reveling the presence of a possible bipolar 
mole cular outflow in NH3 and CO (Bachiller & Cernicharo 
1990 ), having also a cluster of IR sources (Tapia et al. 
19961) . The second region is NGC 6334 1(N)w, which also 
presents a bipolar molecular outflow in SiO (Megeath & 
Ticftrunk |1999| ), is the most intense ammonia source in 
the sky, but is not spatially resolved in our observations. 
It is formed by clumps with diameters of 0.007 pc, masses 
of 0.2 Mq, kinetic temperatures of 24 K and with a total 
mass of about 7600 Mq . This result suggests lar ger m asses 
for this region than estimated by Kuiper et al. ( 1995 ), as- 
suming LTE conditions. It also exihibits velocity gradients 
almost perpendicular to the Galactic plane, and line- width 
gradients parallel to the Galactic plane, showing larger 
line-widths towards the more evolved sources. Finally, the 
region NGC 6334I(N)e is formed by clumps with physical 
conditions similar to those observed in I(N)w, and spread 
through an area larger than the telescope beam. 

The ammonia abundances in both regions are simi- 
lar. Their values indicate that the regions have either ages 
smaller than 10^ years or that depletion effects are impor- 
tant. 
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Among the three regions identified, NGC 63341 is the 
most evolved while NGC 6334 1(N)e is in a very early star 
formation stage. 

Comparison between physical parameters determined 
under LTE and non-LTE conditions showed that the NH3 
column density and ammonia abundance are the only 
quantities which are similar between them. Other param- 
eters, as rotational and kinetic temperatures, exhibit dif- 
ferences that can reach a factor of three. 
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